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Tho rol© of the eubcrltical asaoubly Is gradually being 
enlarged to include studies of transient or dynuniic rosoonsoo. 
Denands for nn inexpensive and quick rothod foi* detorjrsining 
the effect of visiriouo paronoters upon control rosoonao have 
led Investigations into the utilisation of tcio as.:.ewoly for 
such pixrposos. Porhaos tho most iraportant advantage of the 
mbcritlcal assembly is tho ease with viiich exoerlrnonts nay 
bo porfornod. Extra physical »afe^ai‘da and restrictions are 
at a tainitTiUtn, and no extonsivo shielding is required. Tho 
relatively s-Tiall also of tho suboritical asconbly also tends 
to make it easier and leas exnonsivo to vju'y tVi« desired 
parameters. 

Borrowing from sei'vomo chan ism theory and treating the 
assonbly aa a “black box”# many investigators have aptly 
applied tho concept of tho tran?^for function to doocribo tho 
response to a given forcing function. In practice, phase and 
amplitude relations, when considered as functions of the 
driving frequency, are collectively deaijsnated as tho aassmoly 
transfer function. 

Ttie paronetrlo forcing or driving function is of two 
types, reactivity forcing and soux’ce forcing. In general. 
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iwiiotlTity fc^i'otriif l» . e-jroTont®d by t;;» chan:® neutron 
(i^inslty f’u® to a n.-^all limonic variation in tho reactivity, 
thit i", 0 h-.'T'^onic rharco in th» ronctor geonetry uccon- 
plishod by oscillating ecrve nautron abaorbar. Oourc® 
forcing, in contrast to roactlvlty forcing, in no way .u"f«cts 
or chajigoa the aauor*.bly geomotry. It is woroly nn in;>x)30d 
ViC'iatioij in tho ui'is^^ry source. 

The objective of this thesis is to jireaent tho sub- 
crltloal aaso-rbly trseisfor function td:ilch was detominod by 
t.*e tooliniqu® of aou. co forcing. 
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• LVIf> O'’ Tli. Lin;;u\TU-.;^ 

The history of the suocrltioal pile and Its role in the 
devolopmont of the first self -sustaining chain reaction has 
beoorMB legend (1), but ^fork is continuing (2) to nake tho 
aasfmbly a testing vehicle and not noroly an object of in- 
vestigation for itsolf. 

Tne Initial aorvonechanisn concoct was presented by 
Pigott, Crever, and Owens (9 arid 10) at tho Knolls Atonic 
Power’ Laboratory. Because tit© roots of tho kinetics oq;^uation 
are nogative, the pile transfor l*unction was noted as oeing 
the sasre an that for a minitiurj phase shift network. This 
per?T!itted it to be roprenentod by a Bode diagran, k pile 
kinetic sinulator, whoso voltage I’esponso represented the 
pile ncuti'on flux response, was successfully utilised for 
the teat to doterraino tlio stability of an autonatic control 
system and Its ability to nointaln a desired power lovel, 

Prans ($) dotermined su.istantially tho saro typo trans- 
fer function but used a noro sc-phiatieatod approach. By 
definition, the ti'ar.'.sfer function relates sinusoidal varia- 
tions; hence it can oo handled by the standard tecliniquos of 
aervomochanirn theory. The ratio of the Laplace transforza of 
the output to that of the input was found to give the trzxnsfor 
function wore convoniontly . In both the above situations it 
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va** oar*ol'ully notad tUat tho flinctlonn dfterniinf'd woio those 
of a pile slinulator and described tha pile roaponooa only to 
the a sale degree that the or^uatlons utilised doscriood tho 
conditions In the pile. 

It re-inined for Hai’rer, Boyar, and ronicoff (0), basod on 
their nov; classical experlnont on the nrgonn^ C?-2 reactor, 
to justify the uso of the transfer function to describe a I’e- 
actor as a control loop ©Icnont. Since their oxporiT.cnt it 
has becO'^e standard procodui'e to determine the transfer func- 
tion of critical aa^enblios and reactors for staoility crl- 
toi*ia. This method Ims also boon applied to nockupr. of pro- 
posed reactors. Boland, Smith, and hico (3) found the aero 
pov,’Oi' tronafer i*unction of the hPi.-III ruokup asscmoly of the 
il'* :-I, a fast critical nsserably. The results wore used to 
detornlns the best oscillator rods end rieasuroment techniques 
to be used on tho actual eesenbiy, wheai constructed. 

In all the i’oferences found, tho infomatlon prinat*lly 
pertained to a critical assemoly and to tho solution of tho 
transfer function by reactivity forcing. Other than a cur- 
sory discussion in Oiasstono (6, p. 2>2), no unclassified 
reference vaa found which contained any information con- 
cerning tho transfer fimctlon of a subcrltical assembly as 
such or conoornlng source forcing. 
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llio r«sponoo of tb« auijorltical unsossibly to oourco 
forolnc bo dorlved frosj the tixao-dojwsndont gonoral diffu- 
sion o uation, Poj* sicpitcifcy all neutrons axo consldoied 
thermal, end hence the equation can bo wi'itten oa (7, p. 
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v^jore AT 1» the mean velocity of the thoi-nal noutrona* In 
gener’al the synbols in tho derivation are In corraon usogo 16, 
pp, 223-2,^3) > >ut the noanin^ are given in the Li vt of 
Syabola for coiwoniont reference. 

Tlie aouroo tona Is made up of tVxr'eo parts. Bk* first 
portion io tho prompt noutrori« result in>'; frota fission. 
According to the i^erai ago theory this n?.»cnpt; nmitron »out»co 
torta nay bo ©xpi’ossed aa ^ 

Oj, e (1 _ a ) i.„ 0 2 

vlioro is tlMJ total fraction of delayed n^utrom; in any 
one generation, fhe necono pox*tion in the dclayod noxitronn 
which result fr<xr. fission fUid which aro oivoti off duj'ing the 
decay of certain fission products. Oonorally, nix groups ar« 
recognised, each gi’xjup having a dlncx-ot© decay constant 
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Luckllr.fe 

i'raotion of avorKgo group of delayed neutrons 

■'recti n of 1^^ o-*^oup of delayed neutrono 

Conoentratlon < f evaras,o fivonp of delayed 
noutron precursors 

Conoentratlon of group of delayed neutron 
procuraoi a 

Thox'mal neutron diffusion coefficient c«* 

Lffootive rnultipllcfttion factor 

Hultiplication factor for an Infinite assembly 

Reutron diffusion length in the Piodorator cm. 

Avex'age neutron lifetime in a finite assembly soo. 

Kean neutron life tine in an Infinite assembly sec. 

Decay constarst of average group of delayed 

neutron pro cursor 8 aoe.*^ 

Decay com t ant of the ijx'oup of delayed 

neutron precursors s©o.“^ 

Humber of neutrons of thomrial energies 



Desonanco escape probability 

ibormai neutron flux noutrons 

cm.^ sec. 

Kate of production of souroo neutrons neutrons 

sec. 

Kacroscopio absorpti. n czHjassootlon cm."^ 

Tima aeo. 

Fermi ago cra.^ 

Tlieimial neutron ssean velocity cn./soc. 

Fro.jueacy of oscillation radiann 
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and fraction [3 P* H()* '-Ith d*^finod as the concen- 
tration of tiio i^h v'l'oup of delayed neutron oviittcrn, tho not 
rate of formation of the procuraor is ,;lvon by 
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Since the rat® of production of the delayed neutrons is 
numorically equal to tho decay I'ate of tho procuraor, thio 
rate eruala X«C* for the six groups undsi coiiSldora- 

tion. Vith the I’emi age of a delayed neutron assiimod equal 
to that of a proraot neutron, the clowing-dovn nonleakag© 
probability will be e~ . llTils io not strictly ti”ae because 
the energy of the delayed neutrons Is genorally loss than the 
average energy of tho prompt neutrons, but the approximation 
Ig considered reasonable (6, p. 227). The delayed neutron 
soux'Ce term nay bo written as 
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Eq. 4 



where -jo Is defined aa roaonanco escape probability, llio 
tlslrd portion is the external sourco tom vhich will bo 
defined as 



= 3^ + o Eq. 0 

wl^iero Sj la the conpiex onplltudo. 

Por ihrther simplification of tho nathenatica, the space 
and tirao vai’iabloa will be conoidorod aep-.-rablo so tliat tho 
partial differential equations nay be reduced to ordinary 
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dlffei*©nt lal equations (6, p. 221) • 

with substitation of t.quationn 2, 4» and 5 Infcb U;uatlon 
1, tlio rosultinx* o:q?ro3sion is 
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Only first modor will bo considex’od f->r th- ainpliflod solu- 
tion, and hence p. 33)* ^’y definition ia 

the ’^ean lifotine of t jo tnermal noutx*ons in an infinite 
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RoUlim, and it equals 



_ ' . With division by 2.^ and 
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vhero k ir the effoctivo nuitinllca tion factor (ind X la the 
avora,'© thermal neutron llfotlriO, both in a finite aasenoly. 
With tfum roarrtnic^’f^nBt and with substitution of k and JL $ 



it is poasibio to write Equation 7 aa 
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With division by ^ , oollocfcion of tho fcoins containing 

and substitution of = nOT , vboro /ir la conaidored conatarit, 

Equation 10 reducao to 
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Eq, 11 



FVcnn tho above dofinitlona 8 and 9 



but since 
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In the sanwj manner 
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With those aubstitations and division by /IT, Equation 11 
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and with ronrranr,on:ont 
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Vifch sourco forcing, k will bo a constant oecau«© tho 
pil® i;eor.ctry is not chcngod. Tim corroaoonciing forcing 
value of noutron donoity trill b« 
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v^«r® Uq is tho ntoady ooraponont and is the conplex at?5pli- 
tude. In an analogous nannoi* 
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By coT-blnation of E-^uationa 15# 19, and 20# L^untion l8 
beoowoa 



"o * "l" " Ut-1) 



(n^+n^» ) 
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Only onp avorngo groun of delayed neutrons will bo con- 
sidered in order to aako tho equation ooro wieldy. This 
introduces son© additional, but snail # error into the results 
(6, p. 251 )* L'auatlon 21 is the conplete fom, but by the 
definition of a ti'ansfor function only Uio sinusoidal portion 
1b of concern. With separation of tho steady and sinusoidal 
coTponents and vfith tho group notation d 2 *opped, the 
sinusoidal portion of Equation 21 oocones 
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Substitution of Equations 19 and 20 into Equation 15 gives 



k (i 



fe 






(n^+nj^e'* '*"*'* ^^ * 






Bq, 23 



il 



c) 









I 



I << i 






f,- 



l 



•r^i*' ir •(‘■“' V * •• ''■*' * 

^ !L2i£>e^- 

■«tV 



tv k. 



1 



4iiil • f1ll«l « »«4 t-A»j «9i 

•^.'•iri M «»<«i -^n. 4iS:^ • 

Mf T«l itm «( lA 

iHJUaNvvt^ •> « 

» C^AAlft *m.' ^ 'r- I 

<M •»«»• f 

*t . «f , ■ , »• «"£:> l-«i^ 5 i^ 



• 11 ^ «• »A| 



# t»; 



4 



»i 



-Vw jlO 1 1 

Tirr-^ ■ I'*** ' 1 



••'t 



•y' 












• f»; 



4 if ^ 










11 



The sinusoidal portion of ^nation ?3 boconos 
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Substitution of ivC|uatlon 25 into Equation 22 gives 
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^ (k - 1) (k - 1) 
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With the tran^sfer function defined as tho I’Rtio of the 
response to tho sinusoidally varying driving force, the 
desired ratio is 



1 - k 



By rearrangc^^ont Lqiiat.ion 26 boco’^ios 
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Tlxo denoninator of tho right side of r,quation 27 tnay be 
simplifiod to give the desired final forra. 
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Figure 1 i» an overall view of the jrabcritical aoseinbly 
vlth the oast cover ronovod. i*®actor grade graphito 
cylinders 60 inchos long and 7 Inches in diar-etor wore 
mchlnod to a square fom 6 inches across, leaving rounded 
cornora. Taeso were stacked In ten columns nine rows high. 
Additional cyllmlers were cut to a 5-ir^oh high by 6-inoh wide 
foita for the top five rows. The final oil© also, with the 
machined cylinders stac»cod fourteen rows high and ton oolunns 
wide, measures 60 by 60 by 79 inches. The whole array of 
graphite is surrounded by a sandwich of plywood, 0.010 inch 
cadniun, and masonite, bolted in sootionn and bracod with a 
framework. Tho oast and west covers are z*cmovabl© to gain 
access to tho pile. 

Tho graphite is placed on a wooden case constructed ao 
that there are three aocosslble spaces under the pile. These 
spaces extend tho north-south width of the assembly. Two 
aluminum tanks, also extending the width of the pile, wore 
filled with v/ater and placed in the two outermost spaooa to 
raodoj’ato the neutrons and to j’aduco radiations coming di- 
rectly from tho soui-ccs benoath the pile. The central space 
v^as used to house tho oscillator shown in Figure 2. 
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Figure 1. 



Tho subcritlcal asuonbly 
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Tiio oscillator anacnbly 
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Th# holos botweon tho Rmphito blocks, r.''ott3uriiig 
approxlnatoly 1- inches acror.n, pomit insortlon of the 
fuel olenento In various lattice aj’r.uigeTra nt » . /iluninun 
cylinders appr xinately 62 inches long each containing seven 
0.2-lnch uraniu!*! slugs constitute tho fuel olemcnta. 

?ho east cover was in place diu’ing tho tost runs, and 
five holos vere (irillejd on tho cover center lino at dlstanoea 
of one foot, tvo foot, thx'eo feet, foui’ fret, and five feet 
from the base. Foz’ all furthei’ discussion those holos will 
be designated by ordinals with tho hole on© foot above tho 
ba»e being called the first hole, tho hole two feet above 
the ban© being called the second hole, and ao on. VI th tho 
O.Jj-inch lattice, as shown in Figure 1, tho cover holes wore 
in line with ©rapty holes between the graphite blocirs, per- 
mitting a ncuti’on detector to bo inserted into tho pile. 

Tho appai'-atus ixscd waa af^sonblod from x^eadily available 
natorial azid parts. S:<treno accuracy and rofinersent of 
design wore not atteniptod. Five Pu-Bo aoui»co3 cylindrical 
in shape, each laeasuring 1 inch in diameter and 1 3/8 inches 
high, x^ero placed in the source retainer of the oscillator 
shown in Figures 2 and 3« Kodez’ation was provided by a 
pai’affln jacket formed about the sova*co retainer. The 
jacket was made by rollirj(j l/l6-inch aluminum sheet to an 
outside diameter of 7 3/8 inches and cementing it to l/O-inch 
end plates. This container was filled by rioeuis of a filler 



19 



hole v;ith Z$ pounds of noltinc point white paraffin and 

elowly cooled to minialce tho possibility of shi’inlcago voids. 
Approxinatoly throo inches of paraffin ourroundod tho sources. 
The retainer acted as support for both the non- rotating 
paraffin Jacket and tho boai’irii 5 s on which tho oxtsmal 
cylinder rotated. 

The rotating cylinder was for’^od by j’olling l/O-inoh 
aluniinuia shoot to an outside dla^^oter of 0 inches and cutting 
it to a length of 1^ incVios, Tho end plates, also of 
aluciinu'^, woro fi‘om inch stock, shrink fitted and acrov^od 
to a collar fitted on the boarding outer I'un casing. 'HYiq 
cylinder was ^icrowod onto tiio end plates porraitting dis- 
a83o-;'5bly. Tho sin© pattoiTi; was cut fron a 0.010-lnoh thick 
sheet of cadmium to tho dir.cnnlons shown in Figure 3 and ce- 
mented with an aliunlnum based C(nent to tho insido of tho 
rotating cylinder. 

For fi?rtber discusolon position 1 will bo used to refer 
to the rotating cylinder and pattern position whenever point 
A of the flattened pattern in Figure 3 is top dead contor. 
Position 2 vdll bo used to indicate that point B is top dead 
contor. Based on tho difference of readings taken with tho 
cadmlust pattern in tho two extreme positions, it is ostimatod 
that approximately 20 per cent of the souz'ce neutrons reached 
thomal speeds in the paraffin. It is also estimated fi*om 
tho accepted nacroscopic absorption cross section of cadmium 
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ROTATING CYLINDER 
WITH PATTERN 



SOURCE RETAINER 



NONROTATING 
PARAFFIN 
JACKET 




Figure 3* Oscillator and cadmium pattern detail 
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that 93 per cent of tho thermal neutrons reachin ;5 the pattern 
were absorbed. 

An extension of one of tho bearinc collars was ^^ooved 
to hold the leather driving’ belt. The rotation frequency was 
variad by controlling; tho voltage to a ^-Inch Blaclc and Decker 
oloctria drill mounted on a four-foot wooden extension, A 
speed reducer was Interposed to permit speeds down to approxi- 
mately 5 RPM, The oscillator was dealf-pmed for speeds aa 
as 1000 RPM, but no teats wore carried out in this range. 

For maximum signal the neutron detector, a 3/3-inoh 
diameter, 8-lnch long tube filled with enriched DP^ gas, was 
inserted In the first hole to the pile midpoint, A Radiation 
Instrument Development Laboratories, model 200-1, scaler was 
used for the higli voltage and initial amplification. The 
aiapllfiej^-dlaoriisinator jucipor was disconnected and the 
amplifier output taken directly to a Nuclear Chicago model 
I6l5 survey meter with the high voltage rendered Inoperative, 
The output of the survey meter was fed to a type 112 direct 
coupled Tektronix amplifier and thence to one arapllflor and 
pen of a two pen Brush recording system, Tho instrument 
arrangement is pictured In Figure 4, and all the circuitry 
utilized is shovm by block diagram In Figure 5* 

The time applied variation in source strength was re- 
corded by means of brush contacts which were mounted on tho 
oscillator frame and the rotor and which were wired parallel 
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Pigopa Kloctponio instmToentation 
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to a lOK ronietor aa ahown In f’lguro 5» r«Ri.'tor 

vas in aei’ion with a lOOf., reslstoz', and a 110-volt 60- 
cyclo voltage was inproscod across them. Tho second re- 
corder anpllfioi* and pen were also connected across tho lOK. 
resistor. With tho rotor contacts open, a 10-volt 60-eyclo 
Input signal vas plotted by the pen. Whenever the rotor 
contacts closed, this signal was shortod out. The brush 
contact nounted on tho frano was approximatoly 120 degroeo 
t'ron rotor top dead contor while the rotating contact was 
set about at point A on the rotor. It z^as found that 
oscillator speeds up to about 300 kPI. could be measured in 
this nannor. At faster speeds the length of tine tho wav© 
is shorted is only a portion of on© wave cycle and la very 
difficult to read with the rather slow pen response. 

.voadljiga were taken with tho oscillator steady in posi- 
tions I and 2 to dotormino tho maxinun range. This nimulatod 
an infinitely slow oscillation. The difforonce in the 
readings was called i>^eadings v/or© then taken iflth tho 

oscillator rotating at vaa-ious speeds, and for each fre- 
quency the poak to peak: amplitude, called nj|^*, was compzu'ed 
with the base value 

The phase migle was dotormlned by plotting position 1 
of tho rotor on tho tlrso wave and comparing it with the 
position of peeU: flux for each cycle. It was noticed during 
prolininai'y testa that moving the probe steadily In and out 
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of the pilfl by hand produced a recorded v;ave havlno an 
atnplitudo less tnan that ultimately iecordod if fcho probe 
were per'^nlttod to stop li'i the pile foi to ton aeconda. 

The pen, however, sco’^od to reaoond in phase with bho notion. 
It was therefor© docided to correct the roooi’dod laniplitude 
for Che two second tirco constant of th© sum^cj wo ter and to 
accept tho recorded phase angle as oeing appi*oxisiately 
corroct . 

J'-ecRuso the theoretical circuit analysis apoeared to be 
too longthy, it was decided to obtain an arnplitudo correction 
expo ririon tally. A sine wave signal generatoi* was connected 
to the cui'Voy no ter input and adjusted to give a 2-volt out- 
put. Ttio scaler ai^pilficr was not included in the test cir- 
cuit because of bhe internal clrcultiry modifications that 
would havo been entailed with t.^e inclusion, dine© the 
acaloi* amplifier vfas in essence anplifying spike voltages 
with a five nicrosocond resolution, it was considered to have 
nogligiblo effect on the variation of tho rooordod signal 
oraplitude. Tho signal generator frequency was increased 
until the survey meter registered 20,000 counts per nlnuto 
which wag approximately the ss2n« as tho reading obtained with 
the px*ob« inserted to tho pile midpoint in the first hole. 

It had been noted prior that as the oscillator was rotated 
very slowly, tho s?irvoy meter roadin;; varlod between approxi- 
mately 19,000 counts per minute and 21,000 counts p«x' minute* 
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Yho 8if;n£il ganei’afcor was adjuatod until tho survey ivotor 
read th««e valueo, and tho extrencs were rar'.ied on tho dial, 
Tho sinusoidal voi-lation between thoso cxtx'oncis v/aa ohtainod 
by nttacning a yardstick solidly to the signal gonorator 
dial to poiTJit a pendulum action, a heavy weight wac cianped 
at variou."* dictancoa on the yardstick and aiamg to slmulato 
tho different cycling froquoncies. v;ith tho ponduluci hold in 
tho extroT.o positions tho amplitude dii’forence, as oefore, 
sinulatod an Infinitoly slow cycling spood. The ratio of tho 
ob served pcaV: to peak airpiitude at any ono frequ^uicy to tho 
araplitude dlffci’onco detcinuinod above was taifon as the cor- 
i’ootlon fastoi' for that frequency. Tho arplitude corr-ootion 
CEurvo Is plotted in Piguro 6. Tho rtngo of readingn for tho 
correction cui’vo vias aovoroly liaitcd by the pendulum length. 
Because of the physical arrangement, ponoulun Ir.njtna Ions 
thaz:i inches were not foaaiblo. This was net a sorioua 
handicap since shorter yondul’im iongths corr.'ispondod to 
faatoj* cycling speeds which wore beyond the upper values of 
tho ®:E:pcrijr»ntal data, Tho 3® inch limitation of tho yard- 
stick, caused by mounting on tho dial, corrssponded to a 
frequency of 3*59 radians per second using olmplo pondulunj 
theory. It is to bo noted that to obtain a frequoncy of ono 
radian per second a 32 . 2 -foot pendulum would have been 

needed. The correction was applied by dividing by the 
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AMPLITUDE CORRECTION FACTOR 
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Figure 6. AmplitTUd© ratio correction factor vs. frequency 
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corroapondlng correction raoto:- at tho stun© l*»oqu*ncy to give 
tho ojcr>«rlm®ntal values of amplitudo rutlo, plotted 

In Plgu3’« 7. 

f-o«© toat runs ve»re attempted uainjf tho second and third 
holes but tho data I'ore not readable. Apparently the 
statlftical variations v«r« of the sarj© order of magnitude 
as the flux variations, masking tho doairod signal. 
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DIJGU3SI0W Or* .c.;.yJL?S 



The px'lmary results of the experlrnent aro presented 
graphically in 5'lgui*e8 7 8. The data fron which those 

figures wei’e drawn are tabulated in the Appendix. 

Iho OTOoth curves are plots of Lquation 23 for different 
values of tho r-rean nouti^on lifetime in a finite assembly and 
for different values of k. Since noither ^ nor k have been 
determined for the assembly in question, values of k = 0.7 
and 0.8 and - 0.07 socontls and 0.1 seconds were chosen to 
indicate tho ti^ond of the curves with both Ji and k and to 
Indicate the general curve shape. In the theoretical deriva- 
tion the six groups of delayed neutrons were replaced with an 
averego group to reduce tho amount of computation. V/ith tho 
total fraction of tho delayed neutrons equal to 0.0073* an 

average X oqual to O.OO soconds*"^ v;as computed from 

1. fiiXi 

X* P« 230). With tho above paramotoro 



in the Equation 28, the aiJ^litudo ratio and phase anglo v^er-o 
determined for various frequencies, 

Moasui'emonta of the neutron flux in tho first hole woro 
taken for oscillator speeds from approximtely 0,5 to 5 
radians per second. Doth phase angle and amplitudo z’atlo 
readings woro obtained. I’ll© recorded phase angle is plotted 
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Figure 7* Amplitude ratio vs* frequency 
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in Plgiu'fl 8. Tho recorded a’^jjjlltud© ratios woro corrected 
for tbo tvo oocO!k 3 tiiK® constant of th® sxsrvcy m<»tor by tho 
correction factor indicated in Fljjure 6. The data for tho 
anplitude ratio coiTootion factor, tabulated in tho Appendix, 
plots aa a straight line on ssni-logarlthraic paper. The 
range vaa incroaaod by straight lino extrapolation. This 
result in valid only over a llrolted I'ange since for all fre- 
quencies tho factor cannot be greater than unity or less than 
zero. Uonc® tho true factor must deviate soriowl'iat from a 
straight line as it approaches tho oxtrerao values. At tho 
lower frequencies tho correction factor itself bocomoa small, 
so the above deviation will have little effect on tho cor- 
rected amplitude i*atio. At tho higher freauencios where tho 
correction factor boooi^es large in comparison to tho readings 
the effect is najyiif ied, and extreme car© would need to bo 
tal<en to justify any extrapolated values. Fortunately the 
frcqucrcios obtained in the ox?>oriir.ontal correction factor 
tests wore as high as those used in tVi© aiT 5 )litude ratio tests 
and no extrapolated values in tho high frequency rang© had to 
be used. The to plotted in Figure 6 is actually that of the 
pendulim, but it corresponds directly to that of tho 
oscillatoi* sine© in both cases it ie tho frequency of tho 
driving function. Tho fro<|joncy was solved for a simple 
pendulum {without damping) from tVio fornula to = /g/a 

^ O 

lAiox'o g is the gravitational acceleration in foot per aecond*" 
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FREQUENCY - RADIANS / SECOND 



Figure 8 



Phase angle vs* frequency 
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anti is t lo ponduluiB longth in feot. 

An intorofiting result of tho oxporir/ont was the docreano 
in nedlan flux level with increasing CJ . The value of the 

V 

riedian line of tho sinusoidal flux variation for each testing 
fi'equoncy was plotted in Figure 9 as it appoai'od on the pen 
recorder, with tho recorder paper scale as the arbitrary 
ordinate scale. This value corresponds to in the theo- 
retical dovolopnent, and it is tabulated in tho Appendix for 
each frequency tested. In the split tin(^ of Equation 21 into 
the steady and tho sinusoidal conditions, it is seen that 
has no theoretical relation with to . henco tno oooervation 
was not predicted by the theory. The only i'.'gical explana- 
tion lies in tho oscillator design and tho riethod chosen to 
accomplish source forcing. No testa vsoro run to doternino 
the exact cause since tho actual median flux level would have 
little effect on tho cycling amplitude ratio. 

Experimontal uncertainty result o from many factors, A 
very slow rotation trial run was nado with the neutron probe 
within three inches of the oscillator to check tho variation 
in tho source. Tho very s’^eill diat.vrtions noted in the 
plotted slno wave could possibly have coon due to neutron 
scattering, voids in the paraffin, and amplifier drift. 

Perhaps the greatest som*co of error in tho experiment was 
the signal f luotiiationr. due to statistical vai‘iatlons and 
noise, ‘iheso variations obscured and distorted the poalts of 
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Figure 9. Median flvix level vs* frequency 
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Figure 10, Views of typical data 
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the adding to the difficult/ of doterTlnlng both tha 

anpilSuda roitio and the anile. 7\/o sa*noleg of typical 

data, dononoti'ubing roadlnga obtained with both lii.ih and low 
gain, aro shoon la Fi^juro 10. 

7o conpare tho oxperl’wnfcal and calculated values of 
a<«plitU(lo r-itlo ar.il phase anslo, it is only necessary to 
obaorvo the .ssanea of tho oai’voa. It is .seen that tho 
exporirontal pointn lie reasonably close to the tlioorotical 
curve plotted for values of k 0,8 and - 0.0? seconds. 
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Tho trarssfor funotion dot«rminod cxpei'inentnlly agreoa 
qulto wo 11 with that predicted by the theory. Wtillo it was 
not the purpose of the thesis, the oxporlnental points indi- 
cate a mean neutron lifotlno in the assembly of approximately 
0,07 seconds and an effootiy© multiplication factor of 
approximately 0,8. 

Tlie rotor assembly therrnallzed and absorbed a sufficient 
percentage of noutron.s to permit obtaining acceptable data 
from the first hole. In the second and higher holes the 
sinusoidal variation in the neutron flux was reduced to the 
order of magnitude of tho ctatiatlcal and noise variations 
end could not be read. With a larger or more efficient 
neutron detector, a thicker paraffin Jacket, or possibly a 
thicker cadmium pattern, the ampiitiide of the sinusoidal 
vor'iatlons could bo substantially increased permitting 
correspondingly smallor deviations and statistical errors, 

Tlie median level of tho pile neutron flux decreased 
with Increasing cycling speed. This result is not predicted 
by the theory and is p.obably due to tho design of the 
oscillator. 

Although t '!20 experimental results wore not extremely 
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accurate it la felt the data can bo ropi’oduced with 

sufficient accui’acy and consistency to rubatantiato tho 
thoox'otical proillctionn , 
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Ifc is rory apparent that the posoibllltios in cycling 
a subcriticnl iissonbly have not beon fully exploited. It re- 
mains to inprov© on and onlarge the range of readings to much 
faster and much slower cycling 3i>eods. This will rov^uii'e 
some TsodificQtion in tbo oscillator drive mechanian. 

.oftdings with « higher degro© of certainty at tho greatoz' 
speoda can bo ootainod by redceignlng the basic circuitry to 
have a lov?ei' tlr.© con-ctant and perhaps to irxlude a dis- 
criminator to ^N^duce the stQti'*.ticai vari tion. One© this 
has been accors^llahcd it v;ould be cloairablo to obtain the 
transfer function of the same pile by tho toohnic(uo of ro- 
activity forcing and to eo:':pnro the results of th© two 
methods. 

ill tno nuggestlonn above nortain to tsie assembly 
itself as being tho object of invostlgation* With tho 
background knowledge of t^ie transfer function, tho whole 
aroa of investigating the utiliaation of tho subcritlcal 
ass'»mbly as a facility for testing conpononts and for 
advanced design research will have been opened. 
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APPrtfDIX 
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m 




^.6 



’'••fele I 


. ^nnlitudo I'atlo 


correction 


fttotor 8 




1 

Inches 


1 

feet 


m 1 

u> 

(- g/l) rad 


u> 

lans/seoond 


Correction 

amplitude 


orx ection 
factor 


30 


4- 


12 . $7 


3.59 


0.V5 


0.19 




^.oa 


15.5 


3.9U 


0.75 


0.15 


20 


1.66 


19.35 


4.4 


0.6 


0.12 




1 . 2 $ 


25.7 


5.00 


0.35 


0.07 



Table 2 . 


anp litudo ratio 


data 






Run 


CJ 


ni* 


n, * 


n. 




radians/® ocond 


X 


X 


X 






’^1* 


■'1* 


12 


0.49 


0.35 


0.65 


0.95 


13 


0.87 


0.65 


0.65 


0.92 


14 


1.26 


0.5 


0.5 


0.86 


15 


1.73 


0.4 


0.4 


0.89 


16 


2.25 


0.3 


0.3 


0.82 


17 


2.73 


0.2 


0.2 


0.69 


13 


3.31 


0.15 


0.15 


0.65 


19 


3.73 


0.1 


0.1 


0.59 
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Table 3* Phase angle data 



Rxan 


U) 

radians/scoond 


Phase angle -radians 


1 


o.UJi 


0.07 


2 


0.715 


0.27 


3 


1.15 


O. 4 I 4 . 


h 


1.73 


0.53 


5 


1.85 


0.62 


6 


3.31 


0.8S 



Table 4* ^'ledian fliix level# 


to 


Level 


radian s/second 




0.49 


7.3 


0,37 


7.3 


1.26 


7.25 


1.73 


7.2 


2.78 


7.2 


3.31 


7.15 


3.73 


7.3 




6.9 


6.^ 


7.0 


9.67 


7.0 


12.2 


7.0 


13.9 


6.9 


16.4 


6.6 


17.9 


6.3 
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